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Abstract 

We have measured the dijet angular distribution in y^=1.8 TeV pp coUisions 

using the D0 detector. Order QCD predictions are in good agreement with 
the data. At 95% confidence the data exclude models of quark compositeness 
in which the contact interaction scale is below 2 TeV. 



Dijet final states in pp collisions can be produced through quark-quark, quark-gluon 
and gluon-gluon interactions. The angular distributions produced by these processes as 
predicted by theory are similar. Therefore the dijet angular distribution is insensitive to 
the relative weighting of the individual hard scattering process and thus is insulated from 
uncertainties in the parton distribution functions (pdf's). Quantum chromodynamics (QCD) 
parton-parton scattering processes, which are mainly ^-channel exchanges, produce dijet 
angular distributions peaked at small center of mass scattering angles, while many processes 
containing new physics arc more isotropic. Thus, the dijet angular distribution provides an 
excellent test of QCD and a means of searching for new physics such as quark compositeness. 
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We have measured the dijet angular distribution with greater precision over a larger dijet 
invariant mass range and a larger angular range than previous measurements. Next-to- 
leading order (NLO) QCD predictions are available [|I],0 and comparisons can be made 
between the data and both leading order (LO) and NLO predictions. 

The value of the mass scale, A, characterizes the strength of the quark substructure bind- 
ing interactions and the physical size of the composite states. In the regime where A ^ 
is valid, the quarks appear almost pointlike and any quark substructure coupling can be ap- 
proximated by a four-Fermi interaction. With this approximation, the effective Lagrangian 
for a flavor diagonal definite chirality current is 0,^: H =A{2t\ j lS?'\qii^^qii){(iii^^qii) where 
A = ±1, and H= L, R for left or right handed quarks. While this is not the only possible 
contact interaction, it is the only one for which calculations are currently available. Since 
the sign of A is a priori undetermined, limits for constructive interference {A=—l) and 
destructive interference {A=+l) are presented. Previously published results from CDF |^ 
on dijet angular distributions have been compared to the same model in which all quarks 
are composite, yielding 95% confidence limits A^ > 1.8 TeV and A^ > 1.6 TeV on the 
interaction scales. 

The D0 detector, described in detail elsewhere 0, measures jets using uranium-liquid 
argon sampling calorimeters that provide uniform and hermetic coverage over a large range 
of pseudorapidity {\r]\ < 4). Typical transverse segmentation is 0.1 x 0.1 in r/ x 0, where (p 
is the azimuthal angle. 

The data are from the 93 pb~^ sample recorded during the 1994-1995 Tevatron run. 
Events are selected using a multi-level trigger. The first level requires an inelastic collision 
by demanding the coincidence of two hodoscopes on either side of the interaction region. In 
the second level, jet candidates are selected using an array of 40 calorimeter trigger towers 
0.8 X 1.6 in 1] X (j), covering |?7| < 4. Four different trigger criteria are defined, each requiring 
a single trigger tower above a different transverse energy (-Et) threshold. The final level, an 
online software trigger, selects events with a jet above a preset threshold. The Et thresholds 
at which the triggers are > 98% efficient for the rj coverage used in this analysis are 55, 90, 
120 and 175 GeV. 

Jets are reconstructed using a fixed cone algorithm with radius 71= ■\/Arp~+~A0^ = 0.7. 
Calorimeter towers with Et greater than 1.0 GeV are used as seed towers for jet finding [0. 
Jet Et is defined as the sum of the Et in the towers with 71 < 0.7 from the seed tower and a 
new Et weighted center of the jet is calculated. This process is repeated until the jet center 
is stable. When one jet overlaps another, they are merged into a single jet if they share 
more than 50% of the Et of the lower Et jet. Otherwise, they are spht into two distinct 
jets. 

Jet energy calibration is performed in a multi-step process First, the electromagnetic 
energy scale in the central calorimeter is determined by scaling the energies of electrons from 
Z boson decays so that the corrected Z mass agrees with the value measured at LEP 0. 
Next, the jet response of the central calorimeter is measured using momentum conservation 
in a sample of photon + jet events. After determining the jet response for the central 
calorimeter as a function of jet energy, the relative t] dependent jet response is measured 
using both photon -|- jet and dijet events. One jet (photon) is required to be central and 
the jet response is measured as a function of the rj of the other jet. Jets are also corrected 
for out-of-cone showering losses, underlying event, multiple pp interactions, and effects of 
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302 


20 
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1573 
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524 


13 


8789 


175 


>635 


700 


11 


1074 



uranium noise. 

Quality cuts are required for the two leading Ex jets in each event. These cuts eliminate 
spurious jets that arise from noisy calorimeter cells, cosmic rays, and accelerator losses. The 
efficiencies for these cuts are Et and rj dependent and vary between 90% and 97%. 

The dijet invariant mass and kinematic variables of an event are defined using the two 
highest Et jets. The center of mass scattering angle, 9*, and the longitudinal boost, //boost, 
are expressed in terms of the pseudorapidities of these two jets: r7boost=|(^i+'72) and cos6'*= 
tanh?7*, where ?7*=|(^7i — ^12). The dominance of spin-1 gluon exchange gives the angular 
distribution characteristic of Rutherford scattering: dN/dcos9* oc l/sin^(^*/2). To facil- 
itate a comparison with theory, the angular distributions are measured as a function of 
T^_g2|»7*|_ cos6'*|)/(l— I cos6'*|). This definition translates large values of 6* to small 
values of x (e-g- ^* = 90° x=l)- For Rutherford scattering, dN/dx is independent of x- 

The data are selected using a dijet invariant mass (M) threshold and restricting 



the kinematic cuts in order to limit the jets to regions of uniform acceptance pO 
The dijet invariant mass is calculated assuming massless jets and using the expression: 
= 2£'Ti£'T2(cosh {rji — 772) — cos (0i — 02))- Table | shows the average dijet invariant 
mass, maximum x measured, and the number of events for each of four mass ranges. Both 
|?7i| and \t]2\ are required to be less than 3. To maintain uniform acceptance, we also require 

|'7boost| < 1-5. 

Table || shows the dijet angular distribution, (100/A^)(dA^/dx), in bins of Ax = 1 with 
its statistical error in the four mass bins. The jetrad program is used to determine 
the LO and NLO QCD predictions. The jets at NLO are found using the standard jet 
definition which combines two partons into a single jet if they are both within TZ=0.7 of 
their Et weighted center. We require that two partons also be closer than T^sep x 0.7 with 
7?.sep=1.3 [|T^,0. Figure shows the dijet angular distributions normahzed to unit area 
compared to three different theoretical predictions. The dashed and solid curves show the 
LO and NLO predictions for a single choice of renormalization/factorization scale, ^=Et of 
the leading jet. The dotted curve shows the effects of changing the scale to fi=ET/2 for the 
NLO predictions. The LO predictions are insensitive to the renormalization scale, so only 
one scale is shown. With the large value of Xmax, the effects of higher order QCD become 
apparent. The theoretical predictions are clearly sensitive to the order of the calculation 
and to the renormalization scale. The NLO predictions are seen to be in better agreement 
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FIG. 1. Dijet angular distributions for D0 data (points) compared to jetrad for LO (dashed 
line) and NLO predictions with renormalization/factorization scale //=£'t (solid line). The data 
are also compared to jetrad NLO predictions with //=£'t/2 (dotted line). The errors on the data 
points are statistical only. The band at the bottom represents the correlated it 1 cr systematic 
uncertainty. 



with the data than the LO calculations, especially for large x- 

The dominant source of error on the angular distribution is the uncertainty in the r] 
dependence of the calorimeter energy scale, which is found to be less than 2%. Other 
systematic uncertainties, including rj biases in jet reconstruction, multiple pp interactions, r] 
dependent jet quality cut efficiencies, and effects of calorimeter rj and Et smearing are small. 
Because the data distributions are normalized to unit area, uncertainties in the absolute jet 
energy scale are minimal. All systematic uncertainties added in quadrature are shown as 
a band at the bottom of Fig. 1. The effects of a different pdf were examined by replacing 
the default CTEQ3M ^ with CTEQ2MS The calculated angular distribution is 

insensitive to this change. 

Because the currently available NLO calculations do not implement the effects of both 
QCD and quark substructure, possible effects of quark compositeness are determined using a 
LO simulation [Q. The ratio of the LO predictions with compositeness to the LO predictions 
with no compositeness is used to scale the NLO calculations. Figure ^ shows the dijet angular 
distribution for events with M > 635 GeV/c^ compared to theory for different values of the 
compositeness scale. A"*". The largest dijet invariant mass bin is shown because the effects 
of quark compositeness become more pronounced with increasing dijet invariant mass. 

To obtain a compositeness limit, we constructed the variable R^, the ratio of the number 
of events with x < 4 to the number of events with 4 < x < Xmax- The value x = 4 is 



6 



TABLE II. Dijet angular distribution (100/Ar)(dAr/dx) with statistical uncertainties for the 
four mass bins (GcV/c^). 



260< M <425 425< M <475 475< M <635 M >635 
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• D0 Data 

.... A* = 1.6 TeV 

A* = 1.8 TeV 

... A* = 2.0 TeV 
NLO 



/x=E,/2 




Systematic Error 



MASS{GeV/c') > 635 



FIG. 2. Data compared to theory for different compositeness scales. See text for how compos- 
iteness is calculated for NLO predictions. The errors on the points are statistical and the band 
represents the correlated ± 1 a systematic uncertainty. 



chosen to optimize the sensitivity to quark compositeness. Because the angular distribution 
of jets arising from contact interactions is expected to be more isotropic than that for 
QCD interactions, contact interactions will produce more events at small x than QCD and 
therefore will have a larger value of R^. Table |T| shows the experimental ratio for the 
different mass ranges with their statistical and their systematic uncertainties, which are fully 
correlated in mass. Figure ^ exhibits R^ as a function of M for two different renormalization 
scales along with the theoretical predictions for different compositeness scales. Note that 
the two largest dijet invariant mass bins have a lower Xmax value (Table |), and thus a higher 
value of R^ is expected independent of compositeness assumptions. Also shown in Fig. |^ 
are the values for the four degrees of freedom for different values of the compositeness 
scale. The is defined as X^=Z]j,j ^i^ij^^i' where Si is the difference between data and 
theory in each mass bin i. The covariance matrix, V, is defined as V^i=crj^(stat.)+crf (syst.), 
Vij = (Ji(syst.)(Tj(syst.), for i ^ j. For both renormalization scales, the data prefer a model 
without quark compositeness. The data are better fit with ^=Ej-. 

We employed a Bayesian technique [0 to obtain a compositeness scale limit from our 
data, using a Gaussian likelihood function, P(i?^|^), for R^ as a function of dijet invariant 
mass. The compositeness limit depends on the choice of the prior probability distribution, 
P(^). Motivated by the form of the Lagrangian, P{C) is assumed to be either flat in ^=1/A^ 
or ^=1/A'^. Since the dijet angular distribution at NLO is sensitive to the renormalization 
scale, each renormalization scale is treated as a different theory. The likelihood function 
has the form P(i?^|^) = To determine the 95% confidence limit in A, a limit in 
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TABLE III. Values of R^^^ with statistical and fully correlated systematic uncertainties. 



Mass Range (GeV/c^ ) 




Stat. Error 


Syst. Error 


260-425 


0.191 


0.0077 


0.015 


425-475 


0.202 


0.0136 


0.010 


475-635 


0.342 


0.0085 


0.018 


> 635 


0.506 


0.0324 


0.028 



TABLE IV. The 95% confidence limits for the left-handed contact compositeness scale for 
different models. The prior probability distribution is assumed to be flat in 1/A^. 



Compositeness scale 




Ii=Et/2 


A+ 


2.1 TeV 


2.3 TeV 


A- 


2.2 TeV 


2.4 TeV 




2.0 TeV 


2.2 TeV 



is first calculated by requiring that Q(0 P{R^\i')P{i')di' = 0.95 of Q(oo). The hmit 
in ^ is then transformed back into a limit in A. Table |I^ shows the 95% confidence limits 
for the compositeness scale obtained for different choices of models using a prior probability 
distribution which is flat in 1/A^. 

If we vary the models to include constructive interference (A~), or require only up and 
down quarks to be composite (Aud), the 95% confidence limits for the compositeness scale 
change by 0.1 TeV. If the prior probability distribution is assumed to be fiat in 1/A^, the 95% 
confidence limits are reduced by approximately 6%. These limits are valid for either pure 
left- or right-handed contact interactions. Unlike the earlier measurement 0] using Xmax=5, 
the large x explored here gives greater sensitivity to compositeness terms with constructive 
interference than for destructive interference. 

In conclusion, we have measured the dijet angular distribution with greater precision 
over a larger dijet invariant mass range and a larger angular range than previous measure- 
ments. The data distributions are in good agreement with NLO QCD predictions. The 
compositeness limit depends on the choice of the renormalization/factorization scale, the 
model of compositeness, and the choice of the prior probability function. We have presented 
compositeness limits for models with left-handed contact interference. With 95% confidence, 
the interaction scales A+, A~, and A"^^ all exceed 2.0 TeV. 

We thank R. Harris for the use of his program based on Ref . . We also thank the staffs 
at Fermilab and collaborating institutions for their contributions to this work, and acknowl- 
edge support from the Department of Energy and National Science Foundation (U.S.A.), 
Commissariat a L'Energie Atomique (France), State Committee for Science and Technology 
and Ministry for Atomic Energy (Russia), CNPq (Brazil), Departments of Atomic Energy 
and Science and Education (India), Colciencias (Colombia), CONACyT (Mexico), Ministry 
of Education and KOSEF (Korea), and CONICET and UBACyT (Argentina). 
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FIG. 3. as a function of dijet invariant mass for two different renormalization scales. See text 
for how compositeness is calculated for NLO predictions. The inner error bars are the statistical 
uncertainties and the outer error bars include the statistical and systematic uncertainties added in 
quadrature. The values for the four degrees of freedom are shown for the different values of the 
compositeness scale. The D0 data are plotted at the average mass for each mass range. The NLO 
points are offset in mass to allow the data points to be seen. 
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